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We perform a numerical simulation of quantum turbulence produced by thermal counterflow in 
superfluid 4 He by using the vortex filament model with the full Biot-Savart law. The pioneering work 
of Schwarz has two shortcomings: it neglects the non-local terms of the Biot-Savart integral (known 
as the localized induction approximation, LI A) and it employs an unphysical mixing procedure to 
sustain the statistically steady state of turbulence. For the first time we have succeeded in generating 
the statistically steady state under periodic boundary conditions without using the LIA or the 
mixing procedure. This state exhibits the characteristic relation L = 7 2 v^ s between the line-length 
density L and the counterflow relative velocity v ns and there is quantitative agreement between 
the coefficient 7 and some measured values. The parameter 7 and some anisotropy parameters are 
calculated as functions of temperature and the counterflow relative velocity. The numerical results 
obtained using the full Biot-Savart law are compared with those obtained using the LIA. The LIA 
calculation constructs a layered structure of vortices and does not proceed to a turbulent state but 
rather to another anisotropic vortex state; thus, the LIA is not suitable for simulations of turbulence. 



PACS numbers: 

I. INTRODUCTION 

Quantum turbulence 1,2 , which is the disordered mo- 
tion of a tangle of quantized vortices, has been investi- 
gated since pioneering thermal counterflow experiments 
by Vinen in the late 1950s^. Many experimental, the- 
oretical, and numerical studies have advanced our un- 
derstanding of counterflow turbulence. Many researchers 
have recently focused on quantum turbulence near K. 
One of their main motivations is that quantum turbu- 
lence at K exhibits some similarities with turbulence in 
ordinary fluids 1 . A typical example is Kolmogorov's law, 
which is the most important statistical law of classical 
turbulence and which has been numerically confirmed in 
quantum turbulence toe 4 . In contrast, turbulence at fi- 
nite temperatures (e.g., counterflow turbulence) is a form 
of motion peculiar to two-fluid hydrodynamics^ and has 
no direct analog with turbulence in an ordinary viscous 
fluid. For this reason, turbulence at finite temperatures 
is not currently investigated as much as turbulence at 
K. However, the physics of finite- temperature turbu- 
lence is far from fully understood 6 . Recently, particle 
image velocimetry (PIV) experiments 7 9 , which is a tech- 
nique for visualizing fluid flow using small particles, have 
been performed to investigate counterflow of superfluid 
4 He. In order to interpret these experiments, we need to 
understand the microscopic vortex dynamics in thermal 
counterflow. Our current understanding of counterflow 
turbulence is still deficient. Hence, in this introduction 
we briefly review research of counterflow turbulence and 
reveal some important unresolved problems. 

According to the two-fluid model, superfluid 4 He con- 
sists of an intimate mixture of two fluid components: 
a viscous normal fluid and an inviscid superfluid. The 
density and velocity of the normal fluid are respectively 
denoted by p n and v n , while the superfluid density and 



velocity are respectively denoted by p s and v s . The total 
density, p = p s + p n , is approximately temperature inde- 
pendent, but the relative proportions of the normal fluid 
and the superfluid, p s /p and p n /p, depends strongly on 
the temperature. In this system, any rotational motion 
of a superfluid is sustained only by quantized vortices, 
which have the quantum circulation k = h/m^ where h 
is Planck's constant and 777,4 is the mass of a 4 He atom. 

A thermal counterflow, which is internal convection 
produced by a temperature gradient, is explained by this 
two-fluid model. In a counterflow, entropy and heat are 
carried only by the normal fluid component. Hence, if a 
heat current is applied to the closed end of channel, then 
the normal fluid will flow from warmer areas to cooler 
areas while the superfluid will flow in the opposite direc- 
tion to conserve the total mass. In this way, counterflow 
is induced with a relative velocity between the superfluid 
and the normal fluid of v ns = v n — v s . However, super- 
flow becomes dissipative (superfluid turbulence) above a 
certain critical counterflow velocity. The concept of su- 
perfluid turbulence was introduced by Feynman 10 who 
stated that the turbulent state consists of a disordered 
set of quantized vortices, called a vortex tangle. 

This idea was further developed by Vinen. In order to 
describe amplification of a temperature difference at the 
ends of a capillary retaining thermal counterflow, Gorter 
and Mellink introduced some additional interactions be- 
tween the normal fluid and superfluid (mutual friction)ii. 
Through experimental studies of the second-sound at- 
tenuation, Vinen considered this Gorter-Mellink mutual 
friction in relation to the macroscopic dynamics of the 
vortex tangle 3 . Assuming homogeneous superfluid tur- 
bulence (actually, counterflow turbulence is anisotropic), 
Vinen obtained an equation for the evolution of the vor- 
tex line density (VLD) L(i), which we call Vinen's equa- 
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tion: 



II. EQUATIONS OF MOTION 
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where a and \i are temperature-dependent parameters. 
The first term represents the energy injection from the 
normal fluid to the vortices. The second term denotes 
the energy dissipation of vortices due to reconnection be- 
tween vortices. The first and second terms indicate the 
growth and the degeneration of a vortex tangle, respec- 
tively. Therefore, after the growth period of the VLD, 
the vortex tangle enters a statistically steady state. In 
the steady state, the VLD is obtained by setting dL/dt 
equal to zero, which gives 
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where 7 is a temperature-dependent parameter. This 
relation is able to describe well a large number of obser- 
vations of stationary cases^. 

On the other hand, the nonlinear and nonlocal dynam- 
ics of vortices have long delayed progress in achieving 
a microscopic understanding of quantum turbulence. It 
was Schwarz who made a breakthroug h 12 ! 13 . He inves- 
tigated counterflow turbulence using the vortex filament 
model and dynamic scaling 13 . The observable quanti- 
ties obtained by his calculation agree well with the ex- 
perimental results for the steady state of vortex tangles. 
This study confirmed the idea proposed by Feynman that 
superfluid turbulence consists of a quantized vortex tan- 
gle. However, thermal counterflow turbulence is far from 
being perfectly understood. The numerical simulation 
of Schwarz has serious defects. One is that calculations 
are performed under the localized induction approxima- 
tion (LIA), which neglects interactions between vortices. 
Schwarz reported that as a result the layer structure 
is constructed gradually when periodic boundary condi- 
tions are applied. Of course, this behavior differs from 
experimental observations. In order to remedy this, an 
unphysical, artificial mixing procedure was employed, in 
which half the vortices are randomly selected to be ro- 
tated by 90° around the axis defined by the flow velocity. 
It is only this method that enables the steady state to be 
sustained under periodic boundary conditions. These de- 
fects cause us to conjecture that the LIA is unsuitable due 
the absence of interactions between vortices. To under- 
stand counterflow turbulence properly, simulations have 
to be performed using the full Biot-Savart law without 
using the artificial mixing procedure. 

The contents of this paper are as follows. Section I 
I describes the equations of motion of vortices and the 
numerical calculation method. In Sec. IE, we show nu- 
merical simulations of counterflow turbulence by the full 
Biot-Savart law and some physical parameters such as 
the VLD and anisotropy parameters. Section IV com- 
pares the results obtained using the full Biot-Savart law 
with those obtained using the LIA, and shows that the 
LIA is unsuitable. Section V is devoted to conclusions 
and discussions. 



A quantized vortex is represented by a filament passing 
through a fluid and has a definite direction corresponding 
to its vorticity. Except for the thin core region, a super- 
flow velocity field has a classically well-defined meaning 
and can be described by ideal fluid dynamics. The ve- 
locity produced at a point r by a filament is given by the 
Biot-Savart expression: 
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The filament is represented in parametric form s = 
s(£, t), where si refers to a point on the filament and the 
integration is performed along the filament. Helmholtz's 
theorem for a perfect fluid states that the vortex moves 
with the superfluid velocity at the point. Attempting to 
calculate the velocity v w at a point r = s on the filament 
causes the integral diverge as si — » s. To avoid this, we 
divide the velocity s of the filament at the point s into 
two components: 12 



47T 



2(1+1-) 



1/2 



a/4 



a 



— [ 



(si — s) x dsi 

|S1-S|3 ' 

(4) 

where the prime denotes derivatives with respect to the 
arc length £, ao is a cutoff parameter corresponding to the 
vortex core radius, and Z+ and I- are the lengths of the 
two adjacent line elements connected to point s. The first 
term denotes the localized induction field arising from a 
curved line element acting on itself. The second term 
represents the non-local field obtained by performing the 
Biot-Savart integral along the rest of the filament. 

The LIA, which has been used in several studies^ 2 -"— , 
involves neglecting the second non-local term in Eq. (4). 
The equation for the LIA is often written as 



s = /3s / x s". 
Here the coefficient j3 is defined by 
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where c is a constant of order unity and (7 + /_) 1/2 is re- 
placed by the characteristic radius (R) of curvature of 
the vortex lines. In contrast, calculations without the 
LIA are referred to as full Biot-Savart calculations. 

When counterflow is applied, the applied superfluid 
velocity v s is added to v w , and the total velocity so of 
the vortex filament without dissipation is: 
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At finite temperatures the mutual friction due to the in- 
teraction between the vortex core and the normal-fluid 
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~TjK) a a' ~ 
1.3 0.036 0.014 
1.6 0.098 0.016 
1.9 0.21 0.009 
2.1 0.50 -0.03 

TABLE I: Friction coefficients. 



flow v n is taken into account. The velocity of a point s 
is then given byi^ 

s = s + as' x (v n - s ) - a's' x [s' x (v n - s )], (8) 

where a and a' are the temperature-dependent coeffi- 
cients, and so is calculated from Eq. ([7]). Since a' is 
small compared with a, some authors 1 ^ - — neglect a'. 
In this study, we take into account both a' and a. Ta- 
ble IJ presents the mutual fiction parameters used in this 
studyi^. 

The Gross-Pitaevskii (GP) model might be considered 
more appropriate than the vortex filament model. In con- 
trast with the vortex filament model (which is used by 
several groups including us)^~— , the GP model can rep- 
resent phenomena associated with the vortex core includ- 
ing reconnection, nucleation, and annihilation. However, 
no methods have been established for treating mutual 
friction in the GP model, and huge calculations are nec- 
essary to obtain statistics such as the VLD. In contrast, 
the vortex filament model can incorporate the effect of 
mutual friction phenomenologically using the experimen- 
tally observable parameters a and a' ', and it has a lower 
computational cost for dense vortices. For these reasons, 
the vortex filament model is more suitable for numeri- 
cally simulating counterflow turbulence. 

Needless to say, mutual friction plays an important role 
in counterflow turbulence. Let us assume the LI A with 
Eq. (5) and neglect the term with a' . Then, Eqs. (7) 
and (8) are reduced to 

s = j3s' x s" + v s + as' x (v n - v s - j3s' x s"). (9) 

If mutual friction is absent, the dynamics due to only 
the self-induced velocity conserves the total line length of 
vortices. When mutual friction is present and the coun- 
terflow v ns flows against the local self-induced velocity 
j3s' x s", the mutual friction always shrinks the vortex line 
locally. On the other hand, the relative flow along the 
self-induced velocity yields a critical radius of curvature 
given by 

R c - — . (10) 

When the local radius R at a point on a vortex is smaller 
than R Cl the vortex will shrink locally, while the vortex 
will balloon out when R > R c . Thus, mutual friction 
causes the vortex line length to both grow and decay. 
This dual role of mutual friction sustains the steady state 
of counterflow turbulence, ensuring that a highly curved 



structure whose local radius of curvature is less than R c 
will be smoothed out. 

Some important quantities that are useful for charac- 
terizing the vortex tangle are introduced below 13 . The 
VLD is 

L = hJ c ^ < n > 

where the integral is performed at all vortices in the sam- 
ple volume Q . The anisotropy of the vortex tangle that is 
formed under the counterflow v ns is represented by the 
dimensionless parameters 

7 " = ^ ffr-V-h) 2 ]^ ( 12 ) 
I± = ^ z J[l-(s'-T ± ) 2 \d^, (13) 

Here, ry and represent unit vectors parallel and per- 
pendicular to the v ns direction, respectively. Symmetry 
generally yields the relation Iu/2 + I± = 1. If the vortex 
tangle is isotropic, the averages of these parameters are 
I\\ = I± = 2/3 and 1\ = 0. At the other extreme, if the 
tangle consists entirely of curves lying in planes normal 
to v ns , J|| = 1 and I± = 1/2. 

A numerical study of an incompressible Navier-Stokes 
fluid revealed that the close interaction of two vortices 
leads to their reconnection, chiefly because of the vis- 
cous diffusion of the vorticity^. Koplik and Levine di- 
rectly solved the GP equation to show that two close 
quantized vortices reconnect even in a superfluid 17 . Our 
numerical method for vortex filaments cannot represent 
the reconnection process itself. Hence, we reconnect vor- 
tices that pass within the space resolution A£ with unit 
probability. Every vortex initially consists of a string of 
points at regular intervals of A£. When a point on a 
vortex approaches another point on another vortex more 
closely than the fixed space resolution A£, we join these 
two points and reconnect the vortices. This reconnection 
procedure is standard in the vortex filament model, but a 
different procedure is used in some studies^. We discuss 
this in Sec. V. 

In this study, all calculations are performed under the 
following conditions. The numerical space resolution is 
A£ = 8.0 x 10 -4 cm, and the time resolution is At = 1.0 x 
10 -4 s To integrate the equation of motion given by Eq. 
dHj) with respect to time we used the fourth-order Runge- 
Kutta method. The computing box is 0.1x0. 1x0. 1 cm 3 . 
We usually start with an initial vortex configuration of 
six vortex rings, as shown in Fig. QJa). If the vortex ring 
is smaller than R c in Eq. (10), it always shrinks, finally 
disappearing. Hence in our simulation, we make vortices 
that are shorter than Al = 7 x A£ = 5.6 x 10 -3 cm to 
vanish. This cut-off line length is determined to satisfy 
Al < R C1 since a vortex longer than R c has the possibility 
of expanding and causing the VLD to increase. 
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III. NUMERICAL SIMULATION OF 
COUNTERFLOW TURBULENCE 

In this section, we present numerical simulations of 
counterflow turbulence using the full Biot-Savart law un- 
der periodic boundary conditions. Figure [T] shows a typ- 
ical result with the time evolution of the VLD shown in 
Fig. [2] and the anisotropy parameter shown in Fig. 02 
The initial configuration consists of six vortex rings. 

In the first stage (0 < t < 0.4 s), the critical radius 
R c in Eq. (10) determines the vortex destiny. Vortex 
ring sections in which the radius of curvature exceeds 
R Cl expand in the direction perpendicular to v ns through 
mutual friction, while small vortex rings shrink. Thus, 
vortices evolve and become anisotropic, as shown in Figs. 
[2] and 03 At the end of this stage, large vortices appear 
that are comparable to the system size under periodic 
boundary conditions. These vortices survive with a large 
radius of curvature, and continuously generate small vor- 
tices by reconnections in the subsequent stages so that 
they function as "vortex mills" 18 . 

In the second stage (0.4 < t < 2.0s), vortex tan- 
gles undergo continuous evolution despite the decreas- 
ing anisotropy. As vortex rings expand, reconnections 
between vortices occur frequently. Reconnections gener- 
ate vortices with various curvatures, resulting in them 
shrinking and expanding as discussed in the first stage. 
Local sections with a small radius of curvature formed by 
reconnections have an almost isotropic self-induced veloc- 
ity, which prevents the vortices from lying perpendicular 
to v ns , and reduces In. In addition, as the VLD in- 
creases, vortex expansion becomes slower than in the first 
stage because the reconnection distorts vortices, which 
prevents a vortex from smoothly expanding. 

In the third stage (t > 2.0 s), the statistically steady 
state is realized by the competition between the growth 
and decay of a vortex tangle. The growth mechanism is 
still vortex expansion through mutual friction. The decay 
mechanism either creates vortices with local radii of cur- 
vature smaller than R c or the self-induced velocity is ori- 
ented in the opposite direction to v ns after reconnections, 
as discussed in Sec. II. The increasing VLD causes more 
reconnections so that the decay mechanism becomes ef- 
fective. When the VLD has increased sufficiently, the two 
mechanisms begin to compete so that the vortex tangle 
enters the statistically steady state. The LIA calculation 
cannot realize this competition, as discussed in Sec. I 
V, which shows that vortex interaction is essential for 
creating a steady state. 

The steady state is known to exhibit the characteristic 
relation L = 7 2 v 2 s as discussed in Sec. I. Our steady 
states almost satisfy this relation when v ns and L are 
relatively large, as shown in Fig. EJ In Table |TT] we 
show the parameter 7 as a function of T. Our results 
quantitatively agree with the experimental observations 
of Childers and Tough 6,19 . Additionally, there is a crit- 
ical velocity of turbulence, below which vortices disap- 
pear. This critical velocity has been measured in many 




FIG. 1: (Color online). Development of a vortex tangle by the 
full Biot-Savart calculation in a periodic box with a size of 0.1 
cm. Here, the temperature is T = 1.9 K and the counterflow 
velocity v ns = 0.572 cm/s is along the vertical axis, (a) t = 
s, (b) t = 0.05 s, (c) t = 0.5 s, (d) t = 1.0 s, (e) t = 3.0 s, (f) 
t 4.0 s. 



previous studies ^ 21 ! 22 ; it is given by 

2.5 + 1.44a 



u ns,c 



jd 



(15) 



where d is the channel size of the experimental system 
and a is a constant of order unity. In our simulation, 
the system size may be taken to be the size of the peri- 
odic box. Then, Eq. (15) gives v nSiC ~ 0.1cm, which is 
almost consistent with our numerical results. However, 
the temperature dependence of v nSyC should be discussed. 
Equation (15) states that v nSyC should decrease with T, 
which differs from the behavior in Fig. 3J Our numerical 
results show that v nSiC decreases with T below 1.9 K but 
increases at 2.1 K slightly. This is because the strong mu- 
tual friction makes the vortices so anisotropic that they 
cannot form enough reconnections with other vortices, 
and so become degenerate. 



5 



6000 



S 4000 



2000 



v =0.572 cm/s 




FIG. 2: Vortex line density as a function of time for four 
different counterflow velocities. 
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FIG. 4: The steady state vortex line density L(t) as a function 
of the counterflow velocity v ns . The error bars represent the 
standard deviation. 
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show this anisotropy ratio in Fig. [5(b). The anisotropy 
increases with increasing temperature because the mu- 
tual friction increases. The steady state shows slightly 
higher values of I±/I\\ than those obtained by Schwarz. 
This is why the vortex interaction reduces the anisotropy, 
as discussed in Sec. IV. 

All our numerical results are in reasonable agreement 
with experimental results. This means that our model re- 
alistically simulates counterflow turbulence. The steady 
state was obtained without using an artificial mixing pro- 
cedure. 



IV. VALIDITY OF THE LIA 



FIG. 3: Anisotropy parameter as a function of time. Here, 
T — 1.9 K and counterflow velocity v ns = 0.572 cm/s. 



Figure [5(a) shows the anisotropy as a function of 
v ns and T. The anisotropy is almost independent of 
v ns and is dependent on T, in agreement with experi- 
mental observations^ 3 . The anisotropy ratio I±/I\\ has 
been measured experimentally 23 and estimated by nu- 
merical simulation 13 . An isotropic vortex tangle yields 
— 1- ^ the vortex tangle consists entirely of curves 
lying in planes normal to v ns , then I±/I\\ — 1/2. We 



T (K) jnum(s/cm z ) 


7ex P (s/cm^) 


h 


1.3 53.5 


59 


0.738 


1.6 109.6 


93 


0.771 


1.9 140.1 


133 


0.820 


2.1 157.3 




0.901 



TABLE II: Line density coefficients 7 and the anisotropy pa- 
rameter 7|| . jnum and j e x P denote our numerical results and 
experimental results by Childers and Tough^ 9 , respectively. 



Previous studies of counterflow turbulence with the 
LIA have encountered some serious difficulties. Schwarz 
could not obtain the statistically steady state under pe- 
riodic boundary conditions without using a mixing pro- 
cedure, as discussed in Sec. I. Kondaurova et al. could 
obtain the steady state using a different reconnection pro- 
cedure from us, which is discussed in Sec. V. However, 
their values of 7 were three times larger than experimen- 
tally measured values. These difficulties probably arise 
from the LIA. In this section, we compare two calcula- 
tions, namely the LIA (Figs. [6fa),(c)) and the full Biot- 
Savart law (Figs. [6fb),(d)). We run both calculations 
at the same condition, T = 1.6 K and v ns = 0.367 cm/s. 
The time evolution of L(t) and I\\(t) is shown in Figs. [71 
and[8j 

The dynamics of vortices by the LIA is qualitatively 
similar to that with the full Biot-Savart law in the first 
and second stages described in Sec. EL The vortex tan- 
gle appears to enter the steady state for 4 < t < 8 s. 
However, the VLD is about twice that for the full Biot- 
Savart law, as shown in Fig. [7J This difference is due 
to the absence of interaction between vortices; the inter- 
action in general works strongly immediately before and 
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FIG. 5: (a) Anisotropy parameter I\\ as a function of v ns . The 
error bars represent the standard deviation, (b) Anisotropy 
ratio I±/I\\ as a function of temperature obtained by us 
(dots), Schwarz 13 (triangles) and experimentally 23 (vertical 
bars) . 



after reconnections. In the case of the full Biot-Savart 
calculation, through the mutual friction the interaction 
between vortices tends to separate two parallel vortices 
and bring two antiparallel vortices closer to reconnection. 
Reconnection generally creates sharp cusps^ in the vor- 
tex lines, so that the vortices separate rapidly with a 
large self-induced velocity. However, the LI A calculations 
result in very different reconnection dynamics. Recon- 
nection between parallel vortices occurs frequently in the 
LIA, whereas they occur very little in the full Biot-Savart 
calculation due to the interaction. After reconnection 
between parallel vortices specific to the LIA calculation, 
the newly created vortices cannot separate rapidly be- 
cause the radius of curvature is large, as shown schemat- 
ically in Fig. [9j which tends to get the vortices together. 
Also, the resulting vortices are affected little by the de- 
cay mechanism described in Sec. IE. Consequently, the 
average distance of the LIA is smaller and the VLD is 
larger than those of the full Biot-Savart calculation. 

The expansion by mutual friction tends to straighten 
vortices perpendicularly to v ns , but reconnection sup- 




FIG. 6: (a) Side and (c) top views of the LIA calculation, 
(b) side and (d) top views of the full Biot-Savart calculation. 
All figures are at t = 18.6 s. The system is a (0.2cm) 3 cube. 
Applied normal fluid velocity is v ns = 0.55 cm/s. 
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FIG. 7: (color online). Comparison of the time evolution of 
the vortex line density L(t) for the simulation of Fig. [6] 



presses this tendency by creating a low radius of cur- 
vature. However, as discussed above, in the LIA calcu- 
lations, parallel vortices tend to come together due to 
the absence of interaction; hence, vortices cannot create 
a small radius of curvature, which gradually straighten 
and align vortices. Vortices eventually commence form- 
ing a bundle structure composed of parallel straight vor- 
tices, as shown in Figs. [6^a) and (c), which is bunches 
of straight vortices, not turbulence. In this state, the 
vorticity directions of the vortex bundles differ from one 
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FIG. 8: (color online). Comparison of the time evolution of 
the anisotropy parameter Jy (t) for the simulation of Fig. [6] 



layer to another, as shown in Fig. [TUl 



FIG. 9: Reconnect ion of almost parallel vortices. 




V. CONCLUSIONS 

In this study, we investigated thermal counterflow tur- 
bulence using the vortex filament model. The full Biot- 
Savart law was used in our calculation, unlike previous 
studies by Schwarz^. We obtain the statistically steady 
state without using the mixing procedure that Schwarz 
used to sustain the steady state. Our numerical results 
reveal the characteristic relation L = 7 2 v^ s , which has 
been observed in many experiments, and the parameter 
7 agrees with experimentally measured values. Also, the 
anisotropy parameters are in reasonable agreement with 
experimental observations. 

To investigate the validity of the LI A, we compared the 
LIA with the full Biot-Savart calculation. Kondaurova 
et al. mention that the reason why Schwarz encountered 
difficulties is due to the reconnection procedure; they ob- 
tained the steady state using the LIA and a different 
reconnection procedure from us and Schwarai^. They 
assume that reconnections occur only when the vortices 
are expected to cross each other. Probably the details of 
the reconnection procedure are not relevant to statistical 
quantities such as the VLD and the anisotropy param- 
eters. Kondaurova et al. encountered another difficulty 
in that the value of the parameter 7 is very large. Our 
comparison between the LIA and the full Biot-Savart 
reveal that the difficulties encountered by Schwarz and 
Kondaurova et al. originate from the LIA. The steady 
state of turbulence cannot be sustained in LIA calcula- 
tions because most vortices gradually straighten and lie 
in the planes normal to v ns . Since the full Biot-Savart 
calculation includes the interaction between vortices, our 
findings imply that the interaction between vortices is 
essential for turbulence simulations. 

Paoletti et al. succeeded in visualizing counterflow us- 
ing solid hydrogen particles 8 , and they obtained a bi- 
modal distribution for the particle velocity. Paoletti et 
al. expect that particles dragged by the normal fluid and 
particles trapped in the vortex tangle contribute to this 
bimodality, but this still remains to be confirmed. Fur- 
ther study is required to understand this observation. 



FIG. 10: Vorticity direction of bundles in the layer vortices 
shown in Figs. [6ja) and (c). 

The LIA calculation gives very different vortex prop- 
erties from the full Biot-Savart calculation. We can 
thus conclude that the LIA is unsuitable for simulating 
counterflow turbulence. The interaction between vortices 
plays an important role in making turbulence uniform. It 
probably applies to any kind of quantum turbulence. 
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